On the Chemical Mixing Induced by Internal Gravity Waves

Detailed modeling of stellar evolution requires a better understanding of the (magneto)hydrodynamic processes that mix chemical elements and transport angular momentum. Understanding these processes is crucial if we are to accurately interpret observations of chemical abundance anomalies, surface rotation measurements, and asteroseismic data. Here, we use two-dimensional hydrodynamic simulations of the generation and propagation of internal gravity waves in an intermediate-mass star to measure the chemical mixing induced by these waves. We show that such mixing can generally be treated as a diffusive process. We then show that the local diffusion coefficient does not depend on the local fluid velocity, but rather on the wave amplitude. We then use these findings to provide a simple parameterization for this diffusion, which can be incorporated into stellar evolution codes and tested against observations

Reflection and Ducting of Gravity Waves Inside the Sun

Internal gravity waves excited by overshoot at the bottom of the convection zone can be influenced by rotation and by the strong toroidal magnetic field that is likely to be present in the solar tachocline. Using a simple Cartesian model, we show how waves with a vertical component of propagation can be reflected when traveling through a layer containing a horizontal magnetic field with a strength that varies with depth. This interaction can prevent a portion of the downward-traveling wave energy flux from reaching the deep solar interior. If a highly reflecting magnetized layer is located some distance below the convection zone base, a duct or wave guide can be set up, wherein vertical propagation is restricted by successive reflections at the upper and lower boundaries. The presence of both upward- and downward-traveling disturbances inside the duct leads to the existence of a set of horizontally propagating modes that have significantly enhanced amplitudes. We point out that the helical structure of these waves makes them capable of generating an alpha-effect, and briefly consider the possibility that propagation in a shear of sufficient strength could lead to instability, the result of wave growth due to over-reflection.Comment: 23 pages, 5 figures. Accepted for publication in Solar Physic

Fast and Accurate Coarsening Simulation with an Unconditionally Stable Time Step

We present Cahn-Hilliard and Allen-Cahn numerical integration algorithms that are unconditionally stable and so provide significantly faster accuracy-controlled simulation. Our stability analysis is based on Eyre's theorem and unconditional von Neumann stability analysis, both of which we present. Numerical tests confirm the accuracy of the von Neumann approach, which is straightforward and should be widely applicable in phase-field modeling. We show that accuracy can be controlled with an unbounded time step Delta-t that grows with time t as Delta-t ~ t^alpha. We develop a classification scheme for the step exponent alpha and demonstrate that a class of simple linear algorithms gives alpha=1/3. For this class the speed up relative to a fixed time step grows with the linear size of the system as N/log N, and we estimate conservatively that an 8192^2 lattice can be integrated 300 times faster than with the Euler method.Comment: 14 pages, 6 figure

Scaling in Late Stage Spinodal Decomposition with Quenched Disorder

We study the late stages of spinodal decomposition in a Ginzburg-Landau mean field model with quenched disorder. Random spatial dependence in the coupling constants is introduced to model the quenched disorder. The effect of the disorder on the scaling of the structure factor and on the domain growth is investigated in both the zero temperature limit and at finite temperature. In particular, we find that at zero temperature the domain size, $R(t)$, scales with the amplitude, $A$, of the quenched disorder as $R(t) = A^{-\beta} f(t/A^{-\gamma})$ with $\beta \simeq 1.0$ and $\gamma \simeq 3.0$ in two dimensions. We show that $\beta/\gamma = \alpha$, where $\alpha$ is the Lifshitz-Slyosov exponent. At finite temperature, this simple scaling is not observed and we suggest that the scaling also depends on temperature and $A$. We discuss these results in the context of Monte Carlo and cell dynamical models for phase separation in systems with quenched disorder, and propose that in a Monte Carlo simulation the concentration of impurities, $c$, is related to $A$ by $A \sim c^{1/d}$.Comment: RevTex manuscript 5 pages and 5 figures (obtained upon request via email [email protected]

Survival of contact processes on the hierarchical group

We consider contact processes on the hierarchical group, where sites infect other sites at a rate depending on their hierarchical distance, and sites become healthy with a constant recovery rate. If the infection rates decay too fast as a function of the hierarchical distance, then we show that the critical recovery rate is zero. On the other hand, we derive sufficient conditions on the speed of decay of the infection rates for the process to exhibit a nontrivial phase transition between extinction and survival. For our sufficient conditions, we use a coupling argument that compares contact processes on the hierarchical group with freedom two with contact processes on a renormalized lattice. An interesting novelty in this renormalization argument is the use of a result due to Rogers and Pitman on Markov functionals.Comment: Minor changes compared to previous version. Final version. 30 pages. 1 figur

Distributed-order fractional wave equation on a finite domain. Stress relaxation in a rod

We study waves in a rod of finite length with a viscoelastic constitutive equation of fractional distributed-order type for the special choice of weight functions. Prescribing boundary conditions on displacement, we obtain case corresponding to stress relaxation. In solving system of differential and integro-differential equations we use the Laplace transformation in the time domain

Statistics of Atmospheric Correlations

For a large class of quantum systems the statistical properties of their spectrum show remarkable agreement with random matrix predictions. Recent advances show that the scope of random matrix theory is much wider. In this work, we show that the random matrix approach can be beneficially applied to a completely different classical domain, namely, to the empirical correlation matrices obtained from the analysis of the basic atmospheric parameters that characterise the state of atmosphere. We show that the spectrum of atmospheric correlation matrices satisfy the random matrix prescription. In particular, the eigenmodes of the atmospheric empirical correlation matrices that have physical significance are marked by deviations from the eigenvector distribution.Comment: 8 pages, 9 figs, revtex; To appear in Phys. Rev.

Early stage scaling in phase ordering kinetics

A global analysis of the scaling behaviour of a system with a scalar order parameter quenched to zero temperature is obtained by numerical simulation of the Ginzburg-Landau equation with conserved and non conserved order parameter. A rich structure emerges, characterized by early and asymptotic scaling regimes, separated by a crossover. The interplay among different dynamical behaviours is investigated by varying the parameters of the quench and can be interpreted as due to the competition of different dynamical fixed points.Comment: 21 pages, latex, 7 figures available upon request from [email protected]

Ionization degree of the electron-hole plasma in semiconductor quantum wells

The degree of ionization of a nondegenerate two-dimensional electron-hole plasma is calculated using the modified law of mass action, which takes into account all bound and unbound states in a screened Coulomb potential. Application of the variable phase method to this potential allows us to treat scattering and bound states on the same footing. Inclusion of the scattering states leads to a strong deviation from the standard law of mass action. A qualitative difference between mid- and wide-gap semiconductors is demonstrated. For wide-gap semiconductors at room temperature, when the bare exciton binding energy is of the order of T, the equilibrium consists of an almost equal mixture of correlated electron-hole pairs and uncorrelated free carriers.Comment: 22 pages, 6 figure

Scaling and Crossover in the Large-N Model for Growth Kinetics

The dependence of the scaling properties of the structure factor on space dimensionality, range of interaction, initial and final conditions, presence or absence of a conservation law is analysed in the framework of the large-N model for growth kinetics. The variety of asymptotic behaviours is quite rich, including standard scaling, multiscaling and a mixture of the two. The different scaling properties obtained as the parameters are varied are controlled by a structure of fixed points with their domains of attraction. Crossovers arising from the competition between distinct fixed points are explicitely obtained. Temperature fluctuations below the critical temperature are not found to be irrelevant when the order parameter is conserved. The model is solved by integration of the equation of motion for the structure factor and by a renormalization group approach.Comment: 48 pages with 6 figures available upon request, plain LaTe